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We investigated the apsidal motion for the multi-planet systems. In the sim- 
ulations, we found that the two planets of HD 37124, HD 12661, 47 Uma and 
HD 82943 separately undergo apsidal alignment or antialignment. But the com- 
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panions of GJ 876 and v And are only in apsidal lock about 0°. Moreover, we 
obtained the criteria with Laplace-Lagrange secular theory to discern whether a 
pair of planets for a certain system are in libration or circulation. 
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The number of known extrasolar planetary systems are quickly growing (Butler et al. 
2003, Paper I) in recent years. The analysis and orbital fitting of the curves of radial velocity 
can reveal the fact that the main-sequence stars may host one or more planets. Recently, 
Fischer et al. (2003) (Paper II) reported the properties of ten multiple-planet systems (see 
Table 8 in Paper II) and pointed out that most of the systems with the ratio of orbital periods 
less than 5:1 are characterized by mean motion resonance (MMR), such as GJ 876(2:1) 
and HD 82943(2:1). The other important feature for these systems is the apsidal lock of 
the orbiting companions, which indicates that the relative apsidal longitudes of two orbits 
librate about a constant, such that two planets have common time-averaged rate of apsidal 
precession. A number of researchers studied this dynamical mechanism in the planetary 
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systems. The work of Ji et al. (2003a) confirmed that the inner two companions of 55 Cnc 
(Marcy et al. 2002) are in a 3:1 resonance and still experience the apsidal phase-locking. In 
fact, one may be familiar to the presence of the apsidal lock in our solar system, such as 
the Galilean satellites of Io-Europa where the difference of their apsidal longitudes Aw are 
anti-aligned (Lee & Peale 2002). In the case of the pulsar PSR 1257+12, Wolszczan & Frail 
(1992) pointed out that there exists almost a Acu-libration about 180° for the two planets 
in a 3:2 MMR. Then a question arises: are all discovered multiple-planet systems in apsidal 
lock? If so, how frequently does this mechanism take place for these systems? 

According to Laplace-Lagrange (L-L) linear secular theory (Murray & Dermott 1999, 
MD99), if Aw is a libration point, then 360° — Aw could also be the other equilibrium 
due to the symmetry of planetary configurations (also see Beauge et al. 2003). So the 
second question is that if two planets undergo alignment or antialignment, could the mirror 
geometry occur? In this Letter, our goal is to examine the apsidal motion for HD 37124, 
HD 12661, HD 82943, 47 Uma, GJ 876. 

2. Results 

In the present work, we performed a preliminary study of various systems by using N- 
body codes(Ji, Li & Liu 2002). We take the masses of the parent star and companions(with 
sim = 1) from Table 1 for each system. The time step is adopted to be 1% to 2.5% of the 
orbital period of the innermost planet. Here we simply account for coplanar configurations 
for these systems. We remain the values of the semi-major axes unchanged for all orbits. 
The eccentricities and arguments of periapse are generated in the orbital parameters web 
based on the best-fit orbital solutions given the nominal observation errors (see Table l) 1 . 
And the remaining arguments of nodal longitudes and mean anomalies are randomly made 
between 0° and 360°. Finallyfive two-planet systems were individually produced 100 pairs 
of coplanar orbits for integration from 1 Myr to 10 Myr. 

Throughout the paper, A and w denote the mean and periastron longitudes of a planet 
respectively. The subscripts b and c separately represent the inner and outer companions 2 . 
Consequently, the relative apsidal longitude is Aw = w^ — w c . 



1 Here we study the possible motion near the best-fit orbital solutions. This may reveal some important 
dynamical features for a certain investigated system. 

2 For the sake of convenience, we exchange the order of the subscripts for two planets of HD 82943 such 
that the inner planet can be always described by the symbol b. 
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2.1. HD 37124 

Paper I revealed that two companions orbit HD 37124, but indicated that with only 30 
observations the eccentricity of the outer planet was poorly determined, with any eccentricity 
between 0.3 and 0.8 fitting the observations within measurement errors (Butler 2003, private 
communication). In our simulations, we let e c initially be in the interval [0.20, 0.60] in an 
attempt to find likely dynamical constraints for e c . However, all orbits can survive for 10 6 yr, 
indicating the system is extremely stable for any initial e c no greater than 0.60. Moreover, 
we noticed that e c ranges from to 0.62 in secular evolution for all cases, coinciding with the 
outcomes by Paper I and Kiseleva-Eggleton et al.(2002). They pointed out that the system 
self-destructed for the critical initial value of e c above 0.65. In addition, there are two typical 
secular variations of both eccentricities for stable orbits: (a) Case of large oscillations:^ is up 
to 0.85 and e c is less than 0.62.(b)Case of small modulations, both e& and e c are less than 
0.40 (see Figure 1). The dynamical origin of high eccentricity for Case (a) may result from 
Kozai mechanism (Kozai 1962;Wu & Murray 2003) that implies the coupled relationship of 
the maximum eccentricity and minimum inclination. 

Furthermore, it is noteworthy that 43% 3 and 1% of the stable orbits are separately in 
the periapse alignment and antialignment. To the best of our knowledge, such results are a 
new finding for this system. In Figure 1, the upper and low panels separately exhibit the 
variations of eccentricities and relative apsidal longitudes. The results of L-L theory (dash 
lines) confirm the numerical results (thick lines) in the amplitudes and periods of libration 
for both eccentricities and Azu. 



2.2. HD 12661 

The orbital fits revealed that the HD 12661 system is close to 11:2 MMR (Paper 
II;Gozdziewski 2003a) or 6:1 MMR (Gozdziewski & Maciejewski 2003b). With the best- 
fit data from paper II, our simulations exhibit all the orbits are stable for 1 Myr. We observe 
that 46% and 49% of the orbits are in Aro-libration about 0° and 180° respectively. More- 
over, Gozdziewski(2003a) indicated that the possibility of Aw-circulation is much less likely 
than Aro-libration about 0° or 180° in the numerical investigations. Here the two types of 
the apsidal librations almost share an equal opportunity to be acted as one of the dynamical 
mechanisms to remain the stability of HD 12661, because the two planets with moderate 



3 The percentage of the aligned and antialigned cases may be affected by the uncertainty of e c , but imply 
that two planets of this system likely tend to be aligned. 
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eccentricities are always far from each other with isolated orbits and protected from close 
encounters. Still, the fact of A-07-libration indicates that one planet travels through the 
periapse or apoapse when the other moves to the opposite side at the same time. 

The critical arguments for 11:2 MMR 0^ can be written: 



Here we define 0\ = 03 = 2A& — 1 1 A c + 3wb + Qzu c , 02 = 02 = 2A& — 11A C + 2-cub + 7zu c and 
03 = &i — &2 = Aw. Most systems in a 11:2 MMR are only in behavior of #i-libration. A 
particular case is exhibited in Figure 2, which not only ten resonant arguments all librate 
but #3 librates about 180° for 1 Myr. Using octupole-level secular theory, Lee & Peale 
(2003) found that if the ratio of the maximum orbital angular momenta is approximately 
equal to a critical value,then the A-ro-libration about 0° or 180° could happen with large 
amplitude variations of the eccentricities. Additionally, we further found that the status 
of Aro-libration are relevant to the ratio of the initial eccentricities and the initial relative 
apsidal longitudes of two planets (see §3). 



The two companions of 47 Uma is a close analog to the duo of Jupiter and Saturn in 
our solar system in the ratios of the planet masses and orbital periods. Laughlin et al.(2002) 
found the system is in Acc-libration about 0° with low eccentricities of the companions. 
Further, Gozdziewski (2002) used MEGNO technique to investigate in wide ranges of the 
orbital parameters for the system. Using the updated parameters by Paper II, we still found 
this system can experience not only a Acc-libration about 0°, but about 180°. We note 
that all the experiments can last for 10 Myr. The fractions in alignment, antialignment and 
circulation are respectively 20%, 10% and 70%. As well known, Jupiter and Saturn are near 
5:2 MMR, but the Aw circulates at present-day. As for 47 Uma, we did find one of the 
critical arguments 3A& — 8A C + 2wb + 3w c temporarily librates for hundreds of years but Aw 
simultaneously circulates in the secular evolution. However, we argue that 47 Uma may be 
captured into 8:3 MMR in the past during the process of migration due to the planet-disk 
interaction. 



0fc = 2A;, — 11A C + kzub + (9 — k)w c , k = 0...9. 



(1) 



2.3. 



47 Uma 



2.4. 



HD 82943 and GJ 876 



The celebrated nature for HD 82943 and GJ 876 is that the ratio of the orbital periods 
of two companions is respectively close to 2:1. The critical arguments for 2:1 resonance are 
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6i = X b — 2 A c + -cu b , 9 2 = Xb — 2A C + zu c , and the relative apsidal longitude is 6 3 = B\ — 6 2 . In 
the simulations, we found three types of stable orbits for HD 82943: (I)Only Q\ librates about 
0° (Il)Case of alignment, G 1 « 6 2 w 6 3 « 0° (Ill)Case of antialignment, 0i « 180°, # 2 ~ 0°, 
9 3 m 180°. In addition, we noticed that 2%, 3% and 2% of the cases belong to Classes 
I, II and III respectively. The geometry of Class I was also discovered by Gozdziewski & 
Maciejewski (2001). But Classes II and III are quite a new discovery that exhibits the HD 
82943 system can undertake apsidal lock (see Figure 3). Recently, Lee M. H.(2003, private 
communication) pointed out that the best-fit solution 4 makes the orbits closer to be aligned 
than anti-aligned in the study of the 2:1 resonance. We used a semi- analytical (Ji et al. 
2003b) model to investigate the apsidal motion for this system and found that although the 
two planets have high eccentricities up to 0.54 and 0.41, the system remain stable for 10 Myr 
because of both dynamical mechanisms of the 2:1 MMR and apsidal lock. 

Several best-fit orbital solutions for GJ 876 appeared in literature (Marcy et al. 2001; 
Laughlin & Chambers 2001;Rivera & Lissauer 2001). Here we presented numerical results 
based on Laughlin- Chambers fit. Two families of orbits are found to be stable: Classes I 
(14%) and II (7%). In summary, for above two systems, we stress that the 2:1 MMR always 
takes place for stable cases and the apsidal alignment (Laughlin & Chambers 2001 ;Lee & 
Peale 2002) seems to have more likely chances. 

On the other hand, as the two orbits of GJ 876 (or HD 82943) are not well separated, 
repeated interplanetary close approaches may lead to disruption of the systems:79% for GJ 
876 rapidly become unstable, and 93% for HD 82943. 

2.5. Other systems 

Butler et al.(1999) reported the first triple-planet system of v And, consisting of a 51 
Peg-like planet with a 4.6-day orbit and two librating outer planets (Rivera & Lissauer 
2000;Chiang, Tabachnik & Tremaine 2001;Chiang & Murray 2002). With the updated 
data(Paper II), we again examined several tests for integration of 10 6 yr and discovered 
that two outer components can still remain in the apsidal lock about 0°. As for 55 Cnc, we 
revealed that the inner two companions perform the asymmetric Aro-libration about 250° 
or 110° in the secular evolution (Ji et al. 2003a). 

Nagasawa et al.(2003) showed that the two planets of HD 168443 and HD 74156 are 



4 Thc velocities for HD 82943 have not yet been published, and the best-fit orbital parameters taken from 
http: / /obswww.unige.ch. 
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both in circulation. Moreover, we performed additional experiments for HD 38529, but simply 
found the planets circulate. The systems of HD 168443, HD 74156 and HD 38529 are similar 
(Table 2), because each of them has a more massive outer planet and a close inner companion. 
The reason for circulation will be explained in §3. Nevertheless, if one of three systems can 
host an extra middle planet with new observations, there may exist librating planets similar 
to those of 55 Cnc and v And. 



3. Discussion on apsidal motion of Aw 

However, one should make clear the conditions for libration or circulation of Aw for 
two planets. Several theoretical works (Malhotra 2002;Lee & Peale 2003;Nagasawa et al. 
2003;Zhou & Sun 2003)were concentrated on the problem. Beauge et al.(2003) presented 
analytical and numerical outcomes on the existence and location of stable equilibrium solu- 
tions about the systems in mean motion resonance. Here we present our results based on 
L-L theory. Following MD99, we have: 

dwf, e c 

—— = A 12 cos{w b - w c ) - An (2) 

(It t\y 

= ~~^2i cos(ro 6 - w c ) - A 22 (3) 

dAw A A .eb . e c . . , 4 . 

—— = A 22 - A n + (-A 21 A 12 ) cos Aw (4) 

at e c Cb 

where An, A 12 , A 21 and A 22 are initially determined constants (see MD99) and 
the eccentricities of each planet. Here we rewrite the coefficients in equation (4), such that 
Bi = A 22 — An and B 2 = (eb/e c )A 21 — {e c je h )A\ 2 . For hierarchical planetary systems (HD 
168443, HD 74156 and HD 38529, the ratio of the semi-major axes a = a b /a c <^ 1), we 
have | Si | 3> | B 2 \ (see Table 2), thus the part of is dominant and Aw circulates(see also 
Nagasawa et al.2003). According to (4), the critical values of eb/e c for Aw-libration: 



,e_b _ ±(7-l)~ y/( 7 - I) 2 + 4 7 /? 2 (5) 
e c ± 2(5 

where the plus sign in ± corresponds to the case for Acc7-libration about 0°, and the minus 
sign to that of Aro-libration about 180°, with (3 = —bf) 2 (a) /by 2 (a) (by 2 (a) are Laplace 

coefficients, j = 1,2) and 7 = (l/ v /a)(/i c / fibjy/ (1 + A*c) /(I + /•*&) (l^b,c are respectively the 
ratio of the masses of Companion b (c) and host star). Hence, there are two requirements 
for apsidal libration of a pair of planets in the absence of mean motion commensurabilities. 
Firstly, the ratio of initial eccentricities of two planets should not be far from the critical 
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values; secondly, the starting relative apsidal longitude should be satisfied to approach 0° 
or 180°. As a paradigm of HD 37124 (see Figure 1), the initial ratio of the eccentricities 
(e&/e c ) ^ 0.30 is close to (e b /e c ) + = 0.34, and Ag7 ~ 7° ~ 0°, thus the two planets are 
aligned. We can apply the criteria to other systems. 

4. Summary 

In this work, we investigated the apsidal motion for the multi-planet systems:we find 
that the two planets of HD 37124, HD 12661, 47 Uma and HD 82943 separately undergo 
the apsidal alignment or antialignment. But the companions of GJ 876 and v And are 
only in apsidal lock about 0°, which means that no mirror configurations occur for them. 
Moreover, we again used the L-L theory to discern whether two planets of a system librate 
or not. The status of Acc-libration depends on the ratio of their initial eccentricities, semi- 
major axes and planetary masses and original relative apsidal longitudes. And it seems to 
be a selection effect for the companions in the multi-planet systems that are involved in the 
apsidal libration. However, this should be further testified by additional observations of the 
planetary systems. 

We thank G. Laughlin, R.P. Butler, D.A. Fischer, M.H. Lee for informative discus- 
sions and an anonymous referee for good comments to improve the manuscript. This work 
is financially supported by National Natural Science Foundation of China(Grants 10203005, 
10173006, 10233020) and Foundation of Minor Planets of Purple Mountain Observatory. 
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Table 1. The astrocentric orbital parameters of seven multi-planet systems 1 



Planet 


M s tar 


M sin i 


P 


a 


e 






Comments 






(Mo) 


(M Jup ) 


(days) 


(AU) 




(deg) 


(MMR) 


(Align) 


(Anti) 


GJ 876 b a 


0.32 


1.06 


29.995 


0.1294 


0.314(0.02) 


51.8(10) 


2:1 


Yes 


No 


GJ 876 c a 


0.32 


3.39 


62.092 


0.2108 


0.051(0.02) 


40.0(10) 


2:1 


Yes 


No 


HD 82943 b b 


1.05 


0.88 


221.6 


0.73 


0.54(0.05) 


138(13) 


2:1 


Yes 


Yes 


HD 82943 c b 


1.05 


1.63 


444.6 


1.16 


0.41(0.08) 


96(7) 


2:1 


Yes 


Yes 


HD 37124 b c 


0.91 


0.86 


153.0 


0.54 


0.10(0.06) 


97(40) 




Yes 


Yes 


HD 37124 c c 


0.91 


1.01 


1942.0 


2.95 


0.40(0.20) 


265(120) 




Yes 


Yes 


47 UMa b d 


1.03 


2.86 


1079.2 


2.077 


0.05(0.04) 


124.3(12.7) 


8:3? 


Yes 


Yes 


47 UMa c d 


1.03 


1.09 


2845.0 


3.968 


0.00(0.03) 


170.8(15.3) 


8:3? 


Yes 


Yes 


HD 12661 b d 


1.07 


2.30 


263.6 


0.823 


0.35(0.03) 


293.1(5.0) 


11:2/6:1 


Yes 


Yes 


HD 12661 c d 


1.07 


1.57 


1444.5 


2.557 


0.20(0.04) 


162.4(18.5) 


11:2/6:1 


Yes 


Yes 


Ups And b d 


1.30 


0.64 


4.617 


0.058 


0.019(0.20) 


200.0(20.0) 








Ups And c d 


1.30 


1.79 


241.16 


0.805 


0.26(0.03) 


249.1(9.3) 


16:3? 


Yes 


No 


Ups And d d 


1.30 


3.53 


1276.15 


2.543 


0.25(0.03) 


264.2(17.2) 


16:3? 


Yes 


No 


55 Cnc b° 


1.03 


0.83 


14.65 


0.115 


0.03 


90.25 


3:1 


No 


No 


55 Cnc c° 


1.03 


0.20 


44.27 


0.241 


0.41 


46.39 


3:1 


No 


No 


55 Cnc d c 


1.03 


3.69 


4780.0 


5.461 


0.28 


218.91 









1 Recently, Lee & Peale (2003) pointed out the difference between the astrocentric and Jacobi orbits for a hierarchical 
system. The values are taken from the references: (a) Laughlin et al.(2001) (b) http//obswww.unige.ch (c) Butler et 
al. (2003) (d) Fischer et al.(2003) (e) Fischer et al.(2003) and Ji et al.(2003) 
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Table 2. The Atu-circulation for three hierarchical planetary systems 



System 3 " 


M b /M c 


a b /a c 


ef>/e c 


Bi/B 2 


(g+/g-) h 


HD 38529 


0.06 


0.035 


0.81 


18.3 


87.3 


HD 74156 


<0.20 


<0.072 


1.86 


23.6 


>19.0 


HD 168443 


0.45 


0.1027 


2.65 


7026.8 


7.0 



a The data for the masses, eccentricities and semi-major axes 

of two planets for three systems are taken from Fischer et 
al.(2003). 

b Here <?+, g_ are two eigenfrequencies from L-L secu- 
lar solutions. For a = a b /a c <C 1, we have g+/g~ = 
(l/y/a)(fi c /li b )y/(l + // c )/(l + W>) = 7- And g + / g_ ~ 
(1/ y/a)(n c / lib) S> 1, indicating g + is dominant. 
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Fig. 1. — The alignment case for HD 37124. The upper two panels show the secular evolution 
of the eccentricities of Companions b and c. Notice that ranges from 0.08 to 0.14, while e c 
varies from 0.28 to 0.30 (Case b). The bottom panel displays the relative apsidal longitude 
Azu librates about 0° with the amplitude of ±15°. The thick lines denote the numerical 
results (dots for e c ) and the dash lines are those from L-L theory. Both of them are in 
good agreement.The initial values:a fe = 0.54 AU,e 6 = 0.089,fi 6 = 199.22°,^ = 83.38°. a c = 
2.95AU, e c = 0.299,ft c = 344.54°, oj c = 291.10°. 
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2x10 4 4X10 4 6x10 4 8x10 4 1x10 

Time(yr) 



Fig. 2. — The antialignment case for HD 12661. Here Q\ = 2\ — 11A C + 3ti7 b + Qm c , 
62 = 2Ab — 11 A c + 2m b + 7m c and 63 = Am. The signs of diamond and plus, respectively, 
represent Q\ and 62- And the thick line denotes 63, which librates about 180° with small 
amplitudes for 10 5 yr. Notice that 6\ librates about 0°, while 62 does about 180° for the 
same timescale. In fact, ten resonant arguments <p k (k = 0...9) all librate for the same 
timescale for this case. The scenario can be seen for 1 Myr. 
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Time(yr) 

Fig. 3.— The antialignment case for HD 82943 (Class III). Here 9 1 = X b - 2A C + w b , 
@2 — A& — 2A C + zu c and # 3 = Azu. Note that two critical arguments of Q\ (the sign triangle) 
and 9 2 (the sign diamond), respectively, librate about 0° and 180° with slight amplitudes for 
10 4 yr. The thick line represent 6> 3 , librating about 180° in the same timescale. The results 
can extend to 10 Myr. 



